Nuclear Structure corrections for Parity Non-conservation in Atomic Ba+ and Ra 
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We calculate the binding energy, charge radii in ^^®Ba-^'*"Ba and ^^''Ra-^^^Ra using the relativistic 
mean field theory which includes scalar and vector mesons. We then evaluate the nuclear structure 
corrections to the weak charges for a series of isotopes of these atoms using different parameters 
and estimate their uncertainty in the framework of this model. Our results will have important 
implication for the ongoing and planned parity non-conservation experiments and atomic structure 
calculations on Ba+ and Ra"*". 



Parity non-conservation (PNC) in heavy atoms have 
provided an important confirmation 0, H 13 of the 
SU(2)xU(l) electro- weak sector of the Standard Model 
(SM). By combining the results of precision measure- 
ments and calculations using sophisticated many-body 
methods 0, 0, 0|j it is possible to extract the nuclear 
weak charge and compare with its corresponding value in 
the SM. A discrepancy between these two values could 
reveal the possible existence of new physics beyond the 
SM 0. 

As first pointed out by Bouchiat and Bouchiat , the 
matrix element of the PNC Hamiltonian scales as Z^. 
It is primarily because of this reason that heavy atoms 
are considered to be the best candidates for PNC exper- 
iments. A high precision measurement of PNC in atomic 
cesium Q has reduced significantly the uncertainty (< 
1%) in the determination of the nuclear weak charge, 
Qw, oi the Cs nucleus and the deviation from the SM is 
about la 8] . It would also be desirable to consider other 
systems that have the potential to yield accurate values 
of the nuclear weak charge. Ba+ and Ra+ deserve spe- 
cial mention in this context. The transition of interest 
are 6s^Si/2 Sd^Dg/j for Ba+ and 7s^Si/2 Gd^D^/^ 
for Ra"*". An experiment is underway for Ba+ using the 
techniques of ion trapping and laser cooling and another 
has been proposed for Ra+ 0, 0|. Relativistic many- 
body calculations 0, 0| have also been carried out on 
these two ions. 

The experimental result needs input from atomic struc- 
ture calculations involving the interplay of electromag- 
netic and weak interactions. However, the small but 
non-negligible effects of nuclear size must be addressed 
before an interpretation of PNC data in terms of the 
fundamental electro-weak couplings is possible. Thus nu- 
clear structure could become a crucial factor in the in- 
terpretation of PNC experiments of increasing accuracy 
[13L I14L I15L Ho . Il7l | . An extensive discussion on the sen- 
sitivity of atomic PNC and electric dipole moments to 
pos sible new physics has been recently reported in Ref. 

There have been earlier studies to determine nuclear 
structure effects in PNC in atomic cesium using non- 
relativistic potentials [Tsl IT^ as well as relativistic mod- 
els . In this letter we present a relativistic calculation 
of these effects for the Ba and Ra isotopes using the rela- 
tivistic mean field theory (RMF). It is motivated by the 



current efforts to observe PNC in Ba+ and Ra+. 

The RMF theory first proposed by Teller and co- 
workers [Hill 113 and later by Walecka 22] and de- 
veloped by others has been fairly successfully applied to 
both nuclear matter and finite nuclei. The method gives 
good description for binding energies, root mean square 
(rms) radii, quadruple and hexadecapole deformations 
and other nuclear properties not only for the spherical, 
but also for the deformed nuclei. The same parameter set 
of the model also describes well the properties of nuclear 
matter. One of the major attractive features of the RMF 
approach is the incorporation of the spin-orbit interaction 
due to to the presence of the one body Dirac Hamilto- 
nian and the nuclear shell structure automatically arises 
from the nucleon-nucleon interaction via the scalor and 
vector mesons. We can therefore expect the RMF calcu- 
lation to provide useful information on nuclear structure 
corrections to atomic PNC. 

In the Standard Model, the electron-nucleon interac- 
tion is mediated by both the photon and the intermediate 
boson The energy involved in the atomic PNC exper- 
iments are usually only a fraction of an cV, while the mass 
of the Z° is ~ 92 GeV, and so the parity non-conserving 
interaction may be written as a contact interaction. We 
have 
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where B stands for ri (neutron) or p (proton). The first 
term grows coherently with nucleon numbers N and Z. 
The second term together with the anapole moment term 
amounts to at most a few percent of the first term in 
heavy atoms. We shall therefore consider only the first 
term. The effective Hamiltonian becomes 
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where the proton and neutron densities, pp.„(r), are nor- 
malized to unity. We have assumed the Standard Model 
nucleon couplings 
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We need the spatial variation of the electron part iplj^tpe 
over the nucleus, its normalization and its dependence 
on nuclear structure. PNC effects are dominated by s- 
electrons (k = —1) coupled to p-electrons (/t = 1). This 
can be expressed as 



(5) 



where C{Z) contains all atomic-structure effects for a 
point nucleus including many-body correlations, M = 
(0)75V's(0) is the normalization factor for single elec- 
tron and f(r) describes the spatial variation [normalized 
such that /(O) = 1]. It is the integrals 



(6) 



qn,p = / f{r)pn,p{r)d r, 



which determine the effect of the proton and neutron 
distributions on the PNC observables. 

PNC effects are proportional to the matrix element 
between two atomic states i and j, 

O = mpnc\j) = ^C{Z)Af[Qw{N, Z)+Q^-{N,Z)\, 

(7) 

where Qw{N,Z) is the weak charge. For the Standard 
Model, the weak charge takes the form at tree level as 



Qw{N,Z) = -N + Z{l-As,in^ 9w). 



(8) 



The nuclear structure correction Q'^^{N^Z) describes 
the part of the PNC effect that arises from the finite 
nuclear size. In the same approximation as © above 

Q-^^{N,Z) ^ -N{qn-l) + Z{l-^s\x? ew){qp~l)- (9) 

The proton (charge) nuclear form factors needed for qp 
and M are generally well known from the measurements 
of the charge distribution of nuclei close to the stable 
valley and many unstable nuclei as well. The neutron 
nuclear form factor needed for g„ is not well determined 
experimentally and is model dependent. To estimate the 
importance of PNC in nuclear structure, the form factor 
can be approximated to the order of {Za)^ for a sharp 
nuclear surface, and neglecting the electron mass in com- 
parison with the nuclear Coulomb potential Q 7 
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In the above, for a sharp nuclear surface density distri- 
bution, the only relevant parameter is the nuclear radius 
R and (r2") = 3/(2?i + 3)i?2". 

One of the motivations for further improving atomic 
PNC experiments is to test the Standard Model parame- 
ters. After the inclusion of radiative corrections, we begin 
by rewriting Equation (8) and (9) in the form 

QwiN,Z) = 0.9878 X [-iV + Z(l-4.0118x)] 

x(1.0 + 0.00782r),(ll) 



X = 0.23124 ± 0.00017 + 0.0036365 - 0.00258T, (12) 

where x is assumed here to be defined at the mass scale 
Mz by modified minimal subtraction , S' is the param- 
eter characterizing the isospin conserving new quantum 
loop corrections and T characterizing isospin breaking 
corrections. The nuclear structure correction to Qw is 
given by 
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{N,Z) = 0.9878x[-iV(g„-l)-KZ(l-4.0118i)(gp-l)] 

(13) 

The coefficients qn,p defined earlier in equation (6) con- 
tain the nuclear structure effects. We have included the 
intrinsic nucleon structure contributions in evaluating the 
nuclear structure correction. We use [la 
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qp,n= I d^r Pp,„(r) -I- -(r2)/,(p,„)AVp,„/Q)f„ f{r) 

(14) 

where {r'^)i^(p.n) a-re the nucleon weak radii and are 
nucleon weak charges. 

The relativistic Lagrangian density for a nucleon- 
meson many-body system [23, l2^ 
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The field for the cr-meson is denoted by a, that of the lo- 
meson by lo^ and of the isovector p- meson by R^. de- 
notes the electromagnetic field, ipi are the Dirac spinors 
for the nucleons, whose third component of isospin is de- 
noted by T^i. Here g^, g^, gp and e^/47r = 1/137 are 
the coupling constants for a, u, p mesons and photon re- 
spectively. M is the mass of the nucleon and , and 
nip are the masses of the a, lu and p-mesons respectively. 
^f,,.^ j^t^u g^^^ pf^K ^j.g f^^i^ tensors for the UJ^', 
and the photon fields respectively. The field equations for 
mesons and nucleons are obtained from the Lagrangian 
of equation H15|) and can be found in Ref 24] . These are 
nonlinear, coupled partial differential equations, which 
are solved self-consistently. These equations are solved 
by expanding the upper and lower components of the 
Dirac spinors tpi and the boson fields wave functions in 
terms of a deformed harmonic oscillator potential basis. 

The total binding energy of the system is 
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where Epart is the sum of single particle energies of 
the nucleons, E„^ E^^, Ep are the contributions of me- 
son energies, and Ec and Epair are the coulomb and 
pairing energy respectively. We have used the pairing 
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gap defined in Ref. |25j to take pairing in to account. 
Ecm = — |41A~^/^ is the non-relativistic approximation 
for the center-of-mass correction. 

We use the parameter set TMl, NL3 and NL-SH ^ 
for our calculations. We note that in TMl parameter set 
has the non-negative value of the quartic self-coupling 
coefficient gs for the omega mesons. In most of the suc- 
cessful parameter sets the quartic self-coupling term for 
sigma meson is negative, so that the energy spectrum is 
unbounded below. Although in normal cases the solu- 
tions are obtained in local minimum. However, all these 
parameter sets give a good account of various proper- 
ties such as binding energy, compressibility, asymmetric 
energy for nuclear matter. 

We have calculated the binding energies, charge radius 
and shift Sr^ „ and Jr^ „ for the barium and radium iso- 
topes with different parameter sets. The binding energies 
agree in all the cases with the experimental values with 
maximum deviation of 5 to 6 MeV out of a total binding 
energy of 1000 MeV for barium isotopes and similarly 
for radium isotopes the deviation for binding energies 
are around 8 MeV 27]. The charge radii rch agree quite 
well with the fitted values '28'| for the barium isotopes to 
within 1%. The difference between the root mean squre 
radius of neutrons and protons grows almost linearly with 
the neutron number for both the barium and radium iso- 
topes. This difference also depends significantly on the 
theoretical model used. The lack of unambiguous pre- 
cise experimental information on the neutron distribution 
means that one must extrapolate to the desired neutron 
properties. We note that there is essentially no model in- 
dependent experimental information on neutron density 
distributions. We next use these radii to estimate the 
nuclear structure effects in PNC. 

The nuclear structure corrections and the weak charge 
for different isotopes of barium and radium evaluated for 
S = T = for different parameter sets. Here one can 
see that the Qp are constant when the neutron number 
increases. However g„ varies slowly as one increases the 
neutron number. Our RMF calculation gives the nuclear 
structure correction for ^^^Ba, Q"^" = 3.444 for TMl, 
3.453 for NL3 and 3.436 for NL-SH forces. Similarly the 
nuclear structure correction for ^^^Ra are Qyy'^ = 15.240 
for TMl, 15.288 for NL3 and 15.215 for NL-SH parame- 
ters. In Fig I, we have plotted the nuclear structure cor- 
rections versus different isotopes of barium (left panel), 
radium (right panel) for the parameter sets used in our 
calculations. It is seen that the N13 parameter gives a 
higher compared to other parameters for radium 

isotopes. 

We next discuss explicitly the correction to the weak 
charge arising from the difference between the neutron 



and proton distributions. The small difference between 
Qn and Qp has the effect of modifying the effective weak 
charge as 



Qw = Q'^-^"" + AQ'^P (17) 
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FIG. 1: N vs. Qw"^ for barium isotopes(left panel) radium 
isotopes (right panel) in diflferent parameter sets 



where 

AQ'^P = A^(l - qn/qp) (18) 

Assuming the difference by a small parameter, R^/ — 
1 + e, we have 

AQ^V"" - NiZaf {0.221 e)/qp (19) 

Our RMF calculation gives AQ'^^'' ^ 0.252 in TMl, 
AQ'^P = 0.260 in NL3, AQ'^^ = 0.243 in NL-SH 
parameters for "4Ba and AQ^^ = 1.301 in TMl, 
AQ'^P = 1.354 in NL3 and AQ^^^ = 1.274 in NL-SH 
parameters for ^^^Ra. 

In conclusion, we have studied the nuclear weak 
charges for Ba and Ra isotopes. Singly charged ions of 
these atoms have been suggested for possible measure- 
ments of PNC. Our RMF calculation yields AQ'^^/Qw 
of 0.35% for i^-^Ba and 1% for 226Ra. These resuUs will 
have an important bearing on high precision studies of 
PNC in a single isotope or a chain of isotopes of Ba+ and 
Ra+. 
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